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ABSTRACT 

The A star 49 Ceti, at a distance of 61 pc, is unusual in retaining a substantial quantity of molecular 
gas while exhibiting dust properties similar to those of a debris disk. We present resolved observations 
of the disk around 49 Ceti from the Submillimeter Array in the J=2-l rotational transition of CO 
with a resolution of 1.0x1.2 arcsec. The observed emission reveals an extended rotating structure 
viewed approximately edge-on and clear of detectable CO emission out to a distance of ~ 90 AU 
from the star. No 1.3 millimeter continuum emission is detected at a 3cr sensitivity of 2.1 mJy/beam. 
Models of disk structure and chemistry indicate that the inner disk is devoid of molecular gas, while 
the outer gas disk between 40 and 200 AU from the star is dominated by photochemistry from stellar 
and interstellar radiation. We determine parameters for a model that reproduces the basic features 
of the spatially resolved CO J=2-l emission, the spectral energy distribution, and the unresolved CO 
J=3-2 spectrum. We investigate variations in disk chemistry and observable properties for a range of 
structural parameters. 49 Ceti appears to be a rare example of a system in a late stage of transition 
between a gas- rich protoplanetary disk and a tenuous, virtually gas-free debris disk. 
Subject headings: astrochemistry — circumstellar matter — planetary systems: protoplanetary disks 
— starsfindividual (49 Ceti) 



1. INTRODUCTION 

A key to understanding the formation of planetary 
systems lies in characterizing the transitional phase be- 
tween the gas-rich primordial disks found around young 
T Tauri stars and the tenuous, virtually gas-free debris 
disks around their main-sequence counterparts. Unfor- 
tunately, disks in this transitional phase are rare and dif- 
ficult to identify. Dust disks around young stars are com- 
monly identified through the "Vega-excess" phenomenon 
(first o bserved using the In frared Astrono mical Satel- 
lite b y lAumann et alj I1984J ; see review by IZuckerman] 
2001), in which an infared excess over the stellar photo- 
sphere is attributed to reprocessing of optical and ul- 
traviolet starlight by thermally emitting circumstellar 
dust grains. 49 Ceti was fi rst identified in this way by 
ISadakane fc Ni shida (1986). The quantity r = Ljr/L}, i 
is often used as an indicator of the "optical depth" of the 
dust disk, as it provides a rough estimate of the quan- 
tity of optical /u ltraviolet light int ercepted and reemit- 
ted by the dust. Dura et alJl| 19931) correlated the IRAS 
Poi nt Source Catalog with the Yale Bright Star Cata- 
log (jHoffleit fc Jasc hck 1991) and identified three A stars 
with t > 10 -3 , indicative of tenuous, optically thin cir- 
cumstellar dust. Two were the stars 8 Pic and HR4796, 
which are now known to host debris disks. The third 
was 49 Ceti, which unlike the other two defies classi- 
fication as a debris disk because it retains a substan- 
tial quantit y of molecular gas, first observed in the CO 
J=2-l lin e (IZu ckcrman et al]fl995t ) and later confirmed 
in J=3-2 (jDent et al.ll2005f ). At a distance of only 61 pc 
(Hipparcos), it is one of the closest known gas-rich cir- 
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cumstellar disk s, farther only than TW Hydrae (51pc; 
lMamajekl r2005). Its outward similarity to a debris disk, 
combined with the substantial quantity of molecular gas 
still present in the system, suggest that the disk may be 
in an unusual transitional evolutionary phase. 

All three high-r A stars are young: H R 4796A has an 
age of 8 ± 2 Myr (|Stauffer et al.lll995l) a nd 8 Pic has 
been placed at - 20 Myr by iThi et all (|2001bt) . con- 
sistent with the age determinati on of 20 ± 10 Myr by 
iBarrado y Navascues et all (1999). The age of 49 Ceti is 
uncertain due to its isolation; unlike 3 Pic or HR 4796A 
there are no known associated low-mass stars to provide 
a corroborating age estimate. Uura et aH (|1998l ) demon- 
strate that on an HR diagram, all three stars exhibit a 
low luminosity for their color, which is likely attributable 
to their young ages ( ~ 10 Myr). Using the ev olutionary 
tracks of ISiess et ail <j2000f ) . IThi et all (|2001bl ) assign an 
age of 7.8 Myr to 49 Ceti based on its position on the 
HR diagram. 

Few conclusive measurements have been made of the 
dust properties in the 49 Ceti system. HST/NICMOS 
coronographic observations of 49 Ceti failed to de- 
tect any scattered light in the near infrared at r > 
1'.'6 ([Weinberger et al.l Il999f ). Recent subarcsecond- 
scale imaging at m id-infrared wavelengths with Keck 
(jWahhai et al.|[2~007f ) revealed dust emission at 12.5 and 
17.9 /im, extended along a NW-SE axis and apparently 
inclined at an angle of 60°. Simple models of the dust 
emission suggest a radial size segregation of dust grains, 
with a population of very small grains (a ~ 0.1/im) con- 
fined between 30 and 60 AU from the star, and a pop- 
ulation of larger grains (a ~ 15/im) from 60 to 900 AU 
from the star. However, the outer radius of this latter 
component is uncertain due to its dependence on the 
millimeter flux, which is not well determined. There are 
two contradictory single dish measurements of the mil- 
limeter dust emission, both w ith modest signal-to-noise. 
iBockelee-Morvan et all (|1994[ ) report a IRAM 1.2 mm 
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flux of 12.7 ± 2.3 mJy, while ISong et all (|2004) report 
a JCMT/SCUBA 850 /zm flux of 8.2 ± 1.9 mJy. These 
measurements are mutually inconsistent for either a ther- 
mal spectrum (F\ oc A -2 ) or a typical optically thin cir- 
cumstellar disk spectrum (F\ oc A -3 ) in this wavelength 
regime. 

If we accept the lower value of the 850 fim flux and 
make standard assumption s about the dust opacity (e.g. 
iBeckwith fe SargentJ Il99lh . then the total mass of the 
49 Ceti dust disk is 0.1 M 8 . If we compare this to 
other nearby dusty disks at potentially similar stages of 
evolution, we find that 49 Ceti, with an 850 /an flux 
of 8.2 mJy at a distance of 61 pc, has a dust mass 
(oc F^oamd 2 ) approximately 80% that of (3 Pic (104.3 
mJy, 19.3 pc; iHolland et al.lll998f ) but only 0.3% that of 
the typical Herbig A e star HD 169142 (554 mJy, 145 pc; 
ISvlvester et~aTlll996D . Thus the 49 Ceti disk appears to 
have a tenuous dust disk more akin to that of the debris 
disk around (3 Pic than a gas-rich protoplanetary disk. 

Studies of the distribution of gas in the 49 Ceti sys- 
tem have been similarly inconclusive, particularly since 
it is not obvious that a substantial reservoir of molecular 
gas should persist in the strong UV field of an A star 
at this apparently advanced stage. Attempts to detect 
pure rotational transitions of the H2 molecule have re- 
sulted in contradictory reports, with iThi et all ([2001a) 
reporting a marg inal detection using SWS/ISO, which 
IChen et all (120061) did not confirm with Spitzer/IRS ob- 
servations; nor did lCarmona et all (|2007r) detect H 2 emis- 
sion with VLT/CRIRES observations. Models of the 
double-pea k ed JC MT CO J=3-2 line profile observed by 
iDent et ail (|2005[ ) indicated that the gas is likely dis- 
tributed in either a very compact disk with ~ 16° incli- 
nation or a more inclined ring of radius ~ 50 AU and 
inclination ~ 35°. The latter was deemed more consis- 
tent with the dust distribution seen in the mid-infrared, 
although it fails to reproduce the high- velocity wings that 
may be present in the CO J=3-2 line profile. 

In order to obtain spatially resolved information on the 
distribution of material in the system, we observed 49 
Ceti with the Submillimeter Array 4 in the J=2-l transi- 
tion of CO and associated continuum. We detect a rotat- 
ing structure of much greater extent than predicted from 
the single-dish measurements, with a large central region 
devoid of molecular gas emissi on. We also model the 
disk emission using CO STA R (|Kamp fe Bertoldil l2000t 
iKamp fe van Z adclhoff 200lj) , a code that combines thin 
hydrostatic equilibrium models of disks with a rich chem- 
istry network and a detailed heating and cooling balance 
to determine gas properties. The models provide some 
insight into basic properties of the disk, including the re- 
gion of photodissociation of CO in the inner disk and the 
spatial extent of the emission. 

The observations are described in Sj2j and results pre- 
sented in Sj3] In Sj4] we discuss the process undertaken to 
model the data, including the basic model structure, the 
initial conditions for the chemistry, and th e initial model 
adopt ed from the dust emission analysis of Wahha fet al.l 
(2007), as well as adjustments to that fiducial model ne- 

4 The Submillimeter Array is a joint project between the Smith- 
sonian Astrophysical Observatory and the Academia Sinica Insti- 
tute of Astronomy and Astrophysics and is funded by the Smith- 
sonian Institution and the Academia Sinica. 



cessitated by the new observations. The parameter space 
is explored in §4.1[ including an investigation of the vary- 
ing influence of chemistry across the model grid, and ^4.21 
discusses the dust properties in the context of the spec- 
tral energy distribution predicted from the gas model. 
The best-fit model is discussed in §4.31 including an a 
posteriori comparison of the model prediction with the 
observed CO J=3-2 spectrum; inadequacies of the model 
are also noted. The results are discussed in the broader 
context of disk evolution in and a summary is pre- 
sented in Sj6] 

2. OBSERVATIONS 

We observed 49 Ceti with the SMA at 230 GHz dur- 
ing an 11-hour track on the night of October 13, 2006. 
Atmospheric phase was extremely stable, with typical 
phase changes of < 15° between calibrator scans (every 
25 minutes) . Seven antennas were used in the "extended" 
configuration, with projected baselines between 15 and 
130 meters. The primary flux calibrator was Uranus, and 
the passband calibrators were the quasars 3C454.3 and 
J0530+135. Gain calibration was carried out using the 
quasar J0132-169, located just 1.3° from 49 Ceti; the flux 
derived for this quasar was 0.93 Jy. The nearby quasar 
J0006-063 was also included to test the quality of the 
phase transfer from J0132-169. 

Two 2-GHz sidebands separated by 10 GHz were used, 
yielding a continuum sensitivity of 0.7 mJy (la). Spec- 
tral resolution in the line was 0.26 km/s, subsequently 
binned to 2.1 km/s, with rms sensitivity 0.030 Jy in a 
single 2.1 km/s channel. The LSR velocities were con- 
verted to heliocentric using an offset of -9.14 km/s. The 
synthesized naturally weighted beam in the CO J=2-l 
line was l'.'0xl'.'2, at a position angle of -78.6°. Imaging 
was carried out using the MIRIAD software package. 

3. RESULTS AND ANALYSIS 

Figure [T] shows the observed line emission from the re- 
gion around 49 Ceti. Four velocity channels are shown, 
with the velocity indicated by the color of the contour 
lines. The observations are centered on the J2000 coor- 
dinates of 49 Ceti; the star symbol indicates the position 
corrected for the proper motion measured by Hipparcos. 
The maximum signal-to-noise ratio in the line is 8. The 
CO J=2-l emission appears to be in an extended rotat- 
ing structure of > 2" radius, apparently viewed close to 
edge-on. The symmetric distribution of the emission in 
the four velocity channels implies a heliocentric velocity 
near 12.2 km/s, consistent with previous determinations 
of the systemic velocity (10. 5 and 9.9 km/s fo r the disk 
and the star, respectively; see lDent et al.ll2005l and refer- 
ences therein) . No emission is detected outside the range 
of velocities shown. The wide separation of the emission 
peaks, combined with a lack of compact, high-velocity 
emission, suggests that the central regions are clear of 
CO J=2-l emission out to ~90 AU radius (~l'/5), twice 
that of the larger ring predicted from the CO J=3-2 sin- 
gle dish data. Table [T] lists the observed parameters of 
the system. 

Assuming optically thin lines and LTE, the total mass 
in CO probed by the J=2-l transition is given by 

hv-ix A21X2 
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TABLE 1 

Observational parameters for 49 Ceti 



Parameter 



CO(3-2) a T7 CO(2-l) 



13 



CO(2-l) 



continuum 



Rest frequency (GHz) 345.796 230.538 220.399 

Channel width 0.27 km s" 1 2.1 km s -1 8.4 km s -1 

Beam size (FWHM) 14" l'.'0xl'.'2 l'.'0xl'.'2 

PA -78.6° -78.6° 

rms noise (Jy beam" 1 ) 0.22 0.030 0.017 
Dust flux (mjy) 

Peak brightness temperature (K) 0.076±0.008 3.5±0.5 < 0.8 

Integrated intensity (Jy km s" 1 ) 9.5±1.9 2.0±0.3 < 0.5 



230.5 (USB 6 ) 
2x2 GHz 
l'.'0xl'.'2 
-78.6° 
7.0 x 10~ 4 
< 2.1 



49 Ceti 
CO J=2-1 




RA offset (arcsec) 



Fig. 1. — A renzogram of SMA observations of 49 Ceti in the 
CO J=2-l line. The beam size is l'.'0xl'/2, and the position an- 
gle is —79°. Contours are -3, 3, and 5 x 37 mjy/beam (the rms 
noise). The position of 49 Ceti is marked with a star symbol, while 
the gr een line indicates the position angle derived by lWahhai et all 
(t2007T) from mid-IR imaging. The contour colors indicate heliocen- 
tric line-of-sight velocity; the four distinct velocities shown are 9.0, 
11.1, 13.2, and 15.3 km/s, in the order of bluest to reddest channel. 
No emission was detected outside this velocity range. 

where the subscript 21 refers to the CO(2-l) transition, 
F is the integrated flux in the line, d is the distance 
to the source (61 pc; Hipparcos), m is the mass of the 
CO molecule, v is the rest frequency of the transition, 
h is Planck's constant, and X2 = j^ 2 - where N2 is the 
population in the J=2 rotational level while N to t is the 
total CO population. The CO mass calculated using this 
method is 2.2 x 1(T 4 M®. Using the canonical CO/H 2 
ratio of 10 -4 this yields a molecular hydrogen mass of 
2.2 M®, consistent wi th the value of 6.3 x 10~ 3 Mj up = 
2.0 M ffi calculated bv lZuckerman et all (|1995fh 

No continuum emission was detected at this combina- 
tion of resolution and sensitivity. This indicates one of 
two things: either the continuum flux is concentrated at 
the center of the disk but the total flux is too low to be 
detected, or the total flux may be larger but spread over 
many beams, so that the brightness within each beam is 
below our detection threshold. These observations were 
sensitive enough to detect the higher continuum flux re- 



ported by Bockelee-Mo rvan et alj (|1994f ) if it were con- 
centrated within a f ew synthesized bea ms. However, an 
extrapolation of the ISong et all (|2004l ) value for a typi- 
cal circumstellar dust spectrum predicts a lower flux by 
a factor of 6, which is just below the detection threshold. 
The lack of an SMA continuu m detection at 23 GHz 
is therefore inconclusive: if the ISong et aTl ([2004) value 
is correct, we would not expect to detect even centrally 
concentrated emission, and so we cannot constrain the 
spatial extent of dust emission through the nondetection 
at 230 GHz. 

4. DISK MODELING 

In order to gain insight into the physical processes 
at work in the 49 Ceti system, we carried out model- 
ing of the disk with CO STA R (|Kamp fc Bertoldil l205ot 
iKamp fc van Zadelhoffl [200lT) . a code which solves the 
chemical equilibrium simultaneously with a detailed 
heating and cooling balance to determine gas properties 
of circumstellar disks. In the following, the salient fea- 
tures of these models are summarized. The chemistry is 
modeled using a network of 48 different species covering 
the elements H, He, C, O, S, Mg, Si, and Fe. The ele- 
mental abundances and key parameters of these models, 
including the stellar mass, radius, effective temperature, 
surface gravity, and ultraviolet flux, are summarized in 
Table [2j The 48 species are connected through 281 re- 
actions, including cosmic ray chemistry, photochemistry 
and the chemistry of excited H 2 . We compute equilib- 
rium chemistry using a modified Newton-Raphson algo- 
rithm. The solution then only depends on the element 
abundances and not on initial conditions. 

We use the results of dust modeling bv lWahhai et al.1 
(l200l and assume large 30 (im black body grains with 
radiative efficiencies of Q\ = 2ira/\ for A > 2ira and 
Qx = 1 otherwise. These grains are efficient absorbers 
and inefficient emitters, thus achieving dust radiative 
equilibrium temperatures of 

T dust = 324 ( aAim )-°- 2 (r AU )- - 4 K . (2) 

Here, L* and Lq are the stellar and solar luminosity 
respectively, a Mm is the grain size in micron and tau the 
distance from the star in astronomical units. The gas 
temperature is derived from a detailed energy balance 
including the most relevant he ating and cooling processes 
(jKamp fc van Zadelhoffll2001h . 

The radiation field consists of both stellar and interstel- 
lar components. The stellar properties are determined 
by a Kurucz model fit to photometric points collected 
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TABLE 2 

Element abundances and parameters used in the disk 

MODELS 



Parameter™ 


Value 




7.5 x 10~ 2 


A c 


1.3 x 1CT 4 


A 


2.9 x icr 4 




4.2 x 10~ 6 


A Si 


8.0 x 10~ 6 


A S 


1.9 x 10~ 6 


A Fe 


4.3 x 10~ 6 


T eff 


10 000 K 


log g 


4.5 


R* 


1.7 Rq 


M» 


2.3 M 


"DV 


4.68 10 -24 cm -2 H - atom -1 



from the literature (IWahhai et alj 


20071 ISvlvester et all 


19961 iBockelee-Morvan et al.lll994; 


Sontr et al.ll2004h: us- 



ing T c ff=10000 K and log q = 4.5, consistent with the 
values quoted by IChen et al . ( 2006), the derived stellar 
luminosity is L* = 26.1Lq and the radius is 1.7 Kq. 
The spectral energy distribution and Kurucz model are 
plotted in Figure [5J i ncluding dereddening according to 
extinction der i ved by ISvlvester et all (|1996l ) and using a 
iCardelli et al.l (|1989f ) extinction law. The solid line in 
the figure denotes the fit to the photometry of a Kurucz 
stellar atmosphere model at the Hipparcos distance of 61 
pc. The dashed line shows the spectral energy distribu- 
tion of the best-fit model of the outer disk as described 
in t j4.3l The interstellar radiation fiel d in the ultra violet 
is assumed to be 1.2 x 10 cm' ~ 2 s" 1 (lHabindll968l) . 

A basic model of the dust disk was constructed accord- 
ing to the B ayesian analysi s of mi d-infrared emission car- 
ried out by IWahhai et al.l (j2007) . Their model consists 
of an inner disk extending from 30 to 60 AU, composed 
primarily of small grains (a ~ 0.1 fim) with a surface 
density of 5 x 10 -8 g/cm 2 , and an outer disk extending 
from 60 to 900 AU composed of larger grains (a ~ 15 
fim) with a surface density of 3 x 10~ 6 g/cm 2 . They 
derive a surface density distribution for the outer disk 
that is constant with radius, yielding a total disk mass 
of 0.35 Mq. From the mid-IR images, they also deter- 
mine a position angle of 125° ± 10° (indicated in Figure 
[T]) and an inclination of 60° ± 15°. We use this model as a 
starting point for the disk structure, since it reflects the 
best available information on the dust density distribu- 
tion. However, since the molecular gas emission provides 
better constraints on some aspects of disk structure, in- 
cluding the vertical density distribution and the surface 
density structure of the outer disk, we introduce refine- 
ments to this initial model where justified, as described 
in £14.0.21 and £14.11 below. For the large grain population, 
our model uses 30 fim grains instead of 15 /im grains, al- 
though the grain size used in this simple model is highly 
degenerate with other disk properties, as discussed in 

021 

To predict gas properties from this dust model, we 
make two primary assumptions: (1) gas and dust are 
well-mixed, (2) the gas:dust mass ratio is constant. We 
initially assume a constant scale height H—2AU, since 
there is n o information on disk scale height from the dust 
model of lWahhai et~aT1 (l2007h : we also begin by retain- 



ing the inner and outer radi i and radially constan t surface 
density structure from the IWahhai et al.l (|2007h model, 
although these assumptions are modified in i j4.0.2l below. 
Throughout the modeling process, we use the canonical 
gas:dust mass ratio of 100 and assume that the disk is 
embedded in interstellar material of density 10 cm -3 to 
avoid model densities dropping to unrealistically low val- 
ues near the boundaries of the numerical grid. 

To compare our models with the SMA data, 
we use the radiative tr ansfer code RATRAN 
(jHogerheiide fc van der Takl I2000D to generate a 
sky-projected image of the CO J=2-l emission predicted 
for the physical model. We then use the MIRIAD task 
uvmodel to sample the image with the combination of 
spatial frequencies and visibility weights appropriate for 
our SMA data. We allow the inclination and position 
angle of the system to vary in order to best match the 
data. 

HD9672 (49 Cet), d=61 pc, L=26.1 L , E(B-V)=0.08 

" 6 | ' ' ' ' ' ' ' Atlas9(10000^4.5j.7R o ) 




logUA] 

Fig. 2. — Spectral energy distribution (d e-re ddened according to 
extinction derived by Sylvester et al. 1996 and Cardclli ct al. 1989 
extinction law) for 49 Ceti using available optical, infrared, and 
submillimetcr photometry. The solid line denotes a Kurucz stellar 
atmosphere model fitted to the photometry using the Hipparcos 
distance of 61 pc. The dot-dashed line s hows the SED for the 
best-fit model of the outer disk see text of i|4.2l for details. 



4.0.1. Inner Disk 

In the inner disk, inside 60 AU, composed primarily 
of small grains, the stellar radiation field raises the dust 
temperature to 1000-2000 K and dissociates most of the 
molecular gas. In this region, the dominant form of car- 
bon is C + , and even hydrogen is predominantly atomic. 
We therefore ignore the inner disk component in subse- 
quent modeling and focus on reproducing the observed 
CO emission with only the outer disk component. 

This lack of molecular gas in the inner di sk is consis- 
tent w ith t he non-detection of w arm H2 by IChen et alJ 
(|2006h and ICarmona et all ()2007f ). and with the lack of 
high-velocity CO emission in Figure [1] The lack of CO 
emission more than 4.3 km/s from the stellar velocity 
is consistent with an absence of CO within a radius of 
~ 90 AU, for gas in Keplerian rotation around a star of 
2.3 Mffl. 
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4.0.2. Outer Disk 

There are three primary features of the observed CO 
emission from the outer disk that we attempted to repro- 
duce with this modeling effort: (1) the separation of the 
emission peaks in the outer channels (~ 3"), (2) the spa- 
tial extent of the CO emission in all channels, and (3) the 
strength of the emission. Reproducing these features of 
the observed CO emission requires several modifications 
to the best-fit IWahhaj et al.1 (|2007f ) model of the outer 
dust disk, including adjustments to the inner and outer 
radii and a departure from the constant surface density 
prescription. 

At first glance, th e inner radius of 60 AU derived by 
IWahhai et al.l (|2007t ) might seem consistent with the lack 
of emission within 90 AU derived from the missing high- 
velocity wings in our data; however, there is a large re- 
gion at the inner edge of the outer disk subject to pho- 
todissociation by stellar radiation which therefore con- 
tributes little to the CO emission. In order to repro- 
duce the separation of the emission peaks, material is 
required interior to this 60 AU radius. We th erefore take 
the uncertainties in the lWahhai et al.1 (|2007| ) dust distri- 
bution into account and allow the inner disk radius to 
vary. However, moving the inner radius closer than ~ 40 
AU to the star results in high-velocity emission that we 
do not observe in the data, while still producing emission 
peaks wider than observed. We therefore set the disk in- 
ner radius at 40 AU, and then adjust the gas densities to 
further reduce the separation of the emission peaks. 

Increasing the total gas mass leads to an elongated 
morphology with an aspect ratio larger than the obser- 
vations, as the optical depth rises throughout the disk. 
To meet the three criteria of (1) enough gas-phase CO 
near the inner disk edge to reproduce the observed peak 
separation, (2) low enough optical depth in the outer 
parts of the disk to keep the emission from becoming 
more elongated than the data (through photodissociation 
by interstellar UV photons), and (3) maintaining an in- 
ner radius large enough to avoid generating high- velocity 
emission that is not present in the data, we must "pile 
up" material at the inner disk edge to enhance shield- 
ing and concentrate emission. We therefore modify the 
initial assumption of constant surface density as derived 
from the infrared analysis, instead adopting an r _e den- 
sity profile. We simultaneously relax the constant scale 
height assumption, introducing a scale height H that in- 
creases linearly with radius r, with proportionality con- 
stant h = H/r. The full 2-D density structure then be- 
comes n(r, z) — r~ e exp (— z 2 /2H 2 ), where the exponent 
e and scale height constant h are varied to obtain the 
best fit to the CO data. 

The power-law surface density profile results in a much 
better match between the model and the observed emis- 
sion peak separation. It also curbs the elongation of the 
emission to some extent, as the vertical column density 
of the outer disk drops and the material far from the 
star becomes subject to dissociation by interstellar ra- 
diation. However, even steep power law indices for the 
surface density profile do not result in a completely pho- 
toevaporated outer disk and consequently produce emis- 
sion that is much more elongated than observed. In a 
next step, we therefore reduce the outer radius from 900 
to 200 AU. While this is at the lower end of the range al- 



lowed bv lWahhai et alj (|2007f) . their derived outer radius 
was based largely on the uncertain millimeter flux mea- 
surement, and the gas geometry is likely a better probe 
of the disk extent. 

4.1. Grid of Disk Models 

After these initial studies of the outer disk, it became 
clear that several model parameters were ill-constrained 
by previously existing data. Specifically, the disk mass is 
constrained only by the weakly-detected and contradic- 
tory millimeter flux measurements; similarly, the density 
power law index e is ill-determined by the infrared ob- 
servations, which are primarily sensitive to inner disk 
emission. The scale height h is also completely uncon- 
strained by the continuum or single-dish measurements, 
neither of which is sensitive to disk structure in the ver- 
tical direct ion. The disk geome try (PA and inclination) 
quoted by W ahhai et al.1 (p007l ) is also subject to large 
uncertainties, due to the irregular shape of the emission 
observed in the infrared. We therefore attempt to bet- 
ter constrain these disk parameters by using our resolved 
CO gas line observations. Gas lines are generally more 
sensitive than dust emission to temperature and density 
gradients, and can thus provide means to break model 
degeneracies. We ran grids of models for the three struc- 
tural parameters (disk mass, density index, scale height) 
and two geometrical parameters (PA, inclination), find- 
ing the best-fit values by calculating and minimizing a 
X 2 value comparing the model to the observed emission 
from the disk. Due to the computational intensity of the 
calculations necessary to determine the chemistry and 
radiative transfer solutions for each model, we ran only 
a sparsely sampled grid of models. In order to ensure 
that the final model reflects all available observational 
constraints, we centered the grid on the fiducial model 
of []4] and adjusted the parameters only as necessary to 
better reproduce the new CO(2-l) observations, moving 
from coarse to fine grids to ensure adequate exploration 
of the parameter space. We use the modeling primarily 
as a demonstration that the basic features of the observed 
CO emission can be reproduced by a simple azimuthally 
symmetric model of disk structure; the "best-fit" model 
should therefore be viewed as representative of an ini- 
tial understanding of the features of the system rather 
than as a conclusive determination of the disk structural 
parameters. 

4.1.1. CO Chemistry Across the Model Grid 

The CO chemistry is dominated by photodissociation 
in a number of UV bands and thus the abundance of CO 
in each model is mostly dependent on the radial and ver- 
tical column densities being able to shield the stellar and 
interstellar UV radiation respectively. In the following 
we briefly discuss some basic characteristics of the model 
grid. 

The surface density in the models is independent of 
the scale height and hence the radial mass distribution 
in each model can be written as M{R) oc R~ e+3 , where 
M{R) denotes the mass inside a radius R. So, as we 
increase the density power law exponent e, the inner re- 
gion of the outer disk harbors a larger fraction of the 
total mass. The densities in this region of the disk be- 
come higher and hence it is easier to obtain the critical 
column densities necessary for UV shielding in the radial 
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direction. On the other hand, a shallower gradient for 
the density distribution translates into higher densities 
in the outer parts of the disk, thus enhancing the verti- 
cal shielding in the outer disk compared to models with 
high e. None of our models is optically thick in the dust 
continuum, so the UV shielding is mainly H2 shielding of 
the CO bands due to their overlap in wavelengths; CO 
self-shielding also plays a role. 

With this basic picture, we can understand the CO 
chemistry displayed in Fig. [3] as a function of disk mass 
(right column) and density gradient e (center column). 
As the total disk mass is increased, CO first starts to 
build up in the radial direction. It can still be dissoci- 
ated by the vertically impinging interstellar UV radiation 
field in the outer regions of the disk (150-200 AU) un- 
til the disk reaches a mass of ~ 17 M ffi , at which point 
it becomes opaque in the CO bands even in the verti- 
cal direction. A shallow density gradient always leads to 
smaller radial column densities at the same reference ra- 
dius, thus pushing the C + /C/CO transition further out 
in radial distance. In our best-fit model of 13 M@, a 
change in e from 2.5 to 1.1 changes the radius for the 
C+/C/CO front from close to 40 AU to 190 AU. 

The scale height h of the models affects only the verti- 
cal density structure in the models. However, since den- 
sity and chemistry are closely intertwined, it can strongly 
impact the overall radial and vertical structure of the CO 
chemistry. From a comparison of the center panel with 
the bottom left panel in Figure El we see that a factor 
2 lower scale height with respect to the best fit model 
(h = 0.02), enhances the CO abundance in the disk sig- 
nificantly, leading to radial and vertical column densities 
that are more than a factor 10 higher with respect to the 
best fit model. The total CO mass increases by a factor 
of 10 as well, with the integrated emission undergoing a 
corresponding dramatic increase. 

Table [3] displays some key results from a subset of grid 
models such as characteristic radial and vertical CO col- 
umn densities, CO masses and total CO J=2-l line emis- 
sion. For all models in the table, the inner radius is fixed 
at 40 AU and the outer radius at 200 AU. 

4.1.2. From Chemistry to Observables 

The predicted CO J=2-l emission for the models in 
Figure [3] is displayed in Figure [H a comparison of these 
figures illustrates the ways in which differences in chemi- 
cal structure are manifested in the observable properties 
of the CO emission. The CO emission is sampled with 
the same spatial frequencies and visibility weights as the 
SMA data and displayed in renzogram form with the 
same velocity structure as in Figure [1] In order to em- 
phasize the relative structural differences between mod- 
els, the contour levels are 15% of the peak flux for each 
model, with the absolute flux indicated by the thickness 
of the contours, and also printed explicitly at the top of 
each panel. 

The decreased shielding in the inner disk caused by re- 
ducing the density gradient e is visible as a lengthening 
of the emission in the central channels and a widening 
of the emission peaks in the outer channels in the low- 
epsilon model (bottom center panel). Increasing e (top 
center panel) leads to enhanced shielding at the disk in- 
ner edge, causing much higher CO fluxes in the outer 
part of the disk and extremely high contrast between the 



inner and outer velocity channels. 

The primary observable consequence of adjusting the 
mass (right panels, top and bottom) is that the increased 
or decreased shielding from extra gas leads to a cor- 
responding increase or decrease in the total CO flux; 
changes to the shape of the emission are minimal, and 
the primary difference between models of different mass 
over the mass range under consideration is simply in the 
relative brightness of the emission. 

Differences in the scale height of the disk similarly 
manifest as differences in the flux scale; however, de- 
creasing the scale height (bottom left panel) also causes 
greater shielding at the inner disk edge, leading to greater 
elongation of the outer velocity channels and causing the 
inner velocity channels to draw together and overlap as 
the CO flux rises throughout the inner areas of the disk. 
An increase in scale height (top left panel) leads to a 
greater area in the front and back of the disk, projected 
along our line of sight, which increases the flux in the 
central channels and leads to a lower contrast between 
the inner and outer channels of the disk. 

4.2. Spectral Energy Distribution 

After converging initially on a model that was able to 
reproduce the observed CO J=2-l emission, we used that 
model to predict the spectral energy distribution. This 
serves as an a posteriori test of the consistency between 
the gas and dust properties in the models and the avail- 
able observables. 

We integrate over the disk volume to obtain the flux 
as a function of wavelength 

F x = (ira 2 /d 2 ) J J 2nr B\(T dust {r, z))n dust (r, z) Q(\) dz dr , 

where d is the distance to the source and ndust is the num- 
ber density of dust grains in cm' 3 . We assume through- 
out a grain density of 2.5 g/cm 3 . 

While the predicted shape of the spectral energy distri- 
bution matches the observations well, the absolute fluxes 
are too high by a factor of ~ 5. Adjusting the tempera- 
ture of the dust grains alters the shape of the SED curve, 
causing it to deviate from the observed shape; we were 
therefore required to increase the gas:dust ratio from 100 
to 500 in order to reproduce the observed photometry. 
This unusually high ratio is likely an artifact of the sim- 
ple assumptions of the model, since little information is 
available about the dust distribution in this system (and 
none at all from our data). For example, the mass of the 
system is likely not all in 30 /jm grains, and a significant 
fraction of the mass may be in larger grains that con- 
tribute little to the infrared emission. Another possibil- 
ity is that the overall gas:dust ratio is consistent with the 
standard value, but that gas and dust are not well-mixed: 
for example, much of the excess emission may arise from 
the inner edge of the disk, which will be directly illumi- 
nated and heated by the stellar radiation. Resolved ob- 
servations of the dust continuum emission would test this 
hypothesis by placing constraints on the spatial distribu- 
tion of the emitting region. Including effects such as this 
would significantly complicate the model presented here, 
as the H2 formation rate would be affected by varying the 
abundance of the dust on which it forms. In general, the 
dust size and gas:dust mass ratio are strongly related by 
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Fig. 3. — Two-dimensional CO abundances in a subset of disk models. The center panel shows the best-fit model (M = 13Me, e = 2.5, 
h = 0.02), while the rows of models above and below show the effects of incrementing and decrementing, respectively, each of the three 
structural parameters that were allowed to vary during the fitting process: h (left column), e (center column), and M^isfc (right column). 
The values of the parameters shown are h = 0.01, 0.03; e = 1.1, 3.5; and M = 9, 17 MU- 



TABLE 3 

Derived quantities from a subset of the 49 Ceti disk models 



(M e ) 


c 


h 


(10 18 cm- 2 ) 


^(co)3 c 

(10 15 cm" 2 ) 


lvl CO 

(io- 4 M e ) 


7 C o(J=2-l) e 
(Jy km s -1 ) 


13 


2.5 


0.020 


2.76 


4.23 


9.66 


2.6 


9 


2.5 


0.020 


0.32 


1.82 


2.46 


1.2 


17 


2.5 


0.020 


13.5 


9.06 


37.2 


6.9 


13 


3.5 


0.020 


15.1 


4.47 


98.0 


11.7 


13 


1.1 


0.020 


0.13 


0.91 


3.74 


2.3 


13 


2.5 


0.010 


42.8 


78.4 


96.6 


14.5 


13 


2.5 


0.030 


0.12 


2.20 


2.97 


1.5 



the total dust surface required to maintain the observed 
quantity of molecular gas; these are in turn dependent 
on the stellar properties determining dust grain temper- 
atures. None of these dust-dependent quantities are well 
constrained by available data. Given the observations 
available and the extremely simplified dust model, which 
not only neglects the size distribution but also the possi- 
bility of a mixture of compositions and opacities, we use 
the simplest assumption of an altered gas:dust ratio in 
order to conduct a consistency check of the temperature 
and density structure of the gas model. 

Decreasing the total dust mass in the model to match 
the SED reduces the grain surface area for H2 formation. 
Thus molecular hydrogen begins to form at larger radii 
and greater depth, compared to the initial model with 
the canonical gas:dust ratio of 100. As a consequence of 
less effective UV shielding, the total CO mass decreases. 
Hence the total mass of the best-fit model has to be in- 
creased slightly to compensate for the lower molecular 



gas fraction. As a secondary effect, the overall gas tem- 
perature of the dust-depleted model also decreases due 
to the diminished photoelectric heating in the disk. The 
corresponding SED predicted for these parameters is in- 
dicated by the dashed line in Figure [2] The mid- infrared 
flux points are underestimated by this mo del because we 
do not include the inner disk component of Wahhai ct al. 
((2003); as our data provide no constraints on the proper- 
ties of the inner disk, we ignore this component and con- 
centrate on the fit to the outer disk. The flux predicted 
by the model SED is consistent with our own continuum 
upper limit reported in Table [TJ 

4.3. Best-Fit Disk Model 

The center panel of Figure 2] shows the best-fit model 
from the grid, with the minor modifications introduced 
by reproducing simultaneously the spectral energy distri- 
bution. The structural and geometric parameters for this 
model are listed in Table [4] The errors given in the table 
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TABLE 4 

Parameters for Best-Fit Disk Model 



Fig. 4. — CO J=2-l emission predicted for the subset of models shown in Figure [3] sampled with the same spatial frequencies and 
visibility weights as the SMA data in Figure [T] The center panel shows the best-fit model, while the rows of models above and below 
show the effects of incrementing and decrementing, respectively, each of the three structural parameters that we allowed to vary during the 
fitting process: h (left column), e (center column), and M^jsJb (right column). The contour levels are displayed in the upper left corner of 
each panel; they are set at 3 and 5 X 15% of the peak flux for each model. The thickness of the contours is proportional to the absolute 
flux: thicker contours indicate that the source is brighter than the data, while thinner contours indicate that it is fainter than the data. 
The contour levels in the center panel are identical to those in Figure [T] Table [4] gives the full list of parameters for the best-fit model. 

This model reproduces the basic features of the CO 
J=2-l emission well, including the strength of the emis- 
sion, the separation of the emission peaks, and the spatial 
extent of the emission. There are still several important 
differences between the model and the data, however, in- 
cluding (1) an inability to reproduce the changes in po- 
sition angle with radius evident in the data (the "wings" 
of emission extending to the southeast and northwest of 
the position angle axis), and (2) the separation of the 
innermost, low- velocity channels. Both of these may be 
indicative of departures from azimuthal symmetry in the 
disk structure, the former possibly indicating a warp in 
the disk and the latter apparently pointing to a deficit of 

, i . , -, . , . , , , , emission along the minor axis of the disk. In none of our 

arc the approximate L-a uncertainty range interpolated . . ° . . , , ,. . . , 

from the % grid models were we able to reproduce the wide separation be- 
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Fig. 5. — Two-dimensional gas temperatures in the best fit 
disk model (M disk = 13 M®, e = 2.5, /i = 0.02. Shown are the 
most important heating (top panel) and cooling (middle panel) 
processes as well as the gas temperature (bottom panel). The dust 
temperature, which depends only on radius, is overlaid in white 
contour lines (steps of 10 K). 

tween the inner channels; while the signal-to-noise ratio 
in these channels is low, the observed CO morphology is 
difficult to reproduce in detail with a simple, azimuthally 
symmetric disk model. The CO emission for this best- 
fit model is optically thin in both the J=2-f and J=3-2 
transitions, even for the edge-on disk orientation, and 
therefore traces the full column density of disk material. 

The densities in the disk are too low for efficient gas- 
dust coupling and thus the gas finds its own equilibrium 
temperature determined mainly by photoelectric heating 
and line cooling. The most important cooling lines from 
surface to midplane are [C n] , [O i] , and H 2 . CO abun- 
dances are only high in a region between 45 and 70 AU 
(Fig- El). Outside that region, CO cooling is less impor- 
tant for the energy balance. Fig. [5] summarizes the most 
important heating and cooling processes and also shows 
the resulting gas temperature structure. The disk surface 
stays well below 100 K due to efficient fine structure line 
cooling. The molecular cooling is however less efficient in 
competing wit h the photoelectric heating fro m the large 
silicate grains (jKamp fc van Zad clhoff 2001), leading to 
temperatures of a few 100 K in the disk interior. 

In order to test the robustness of the best-fit model 
to the gas properties, we used this model to predict the 
CO J=3-2 spectrum. It compar e s fav orably with the 
spectrum observed bv lDent et alj (|2005[ ). reproduced in 
Figure[6] The heavy solid line shows the J=3-2 spectrum 
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Fig. 6. — CO J=3-2 spectrum predicted for the model that pro- 
vides the best fit to the resolved J =2- 1 emission (heavy solid line), 
compared with the IDent et al.l (120051 ") JCMT CO J=3-2 spectrum 
(light solid line). The x-axis shows heliocentric velocity while the 
y-axis gives the JCMT main beam brightness temperature. 

predicted from the best-fit disk model, while the light 
solid line shows the observed JCMT spectrum. Although 
the observed spectrum is noisy and likely subject to an 
absolute calibration uncertainty, the overall agreement 
is within ~ 30%, which is very good given that the CO 
J=3-2 spectrum was not used a priori to determine these 
model parameters. 

5. DISCUSSION 

The processes determining the amount and distribu- 
tion of gas and dust in transition disks like that around 
49 Ceti are the same processes that shape the features 
of emergent planetary systems around these young stars. 
Resolved observations of individual disks in this phase 
are desirable to address such basic questions as when in 
the lifetime of a star its disk disperses, whether the gas 
clears before the dust, and whether the disk clears from 
the center or in a radially invariant manner. 

In the 49 Ceti system, the infrared d ust properties 
appe ar similar to those of a debris disk l|Wahhai et al.l 
2007). Yet observations presented here indicate that a 
substantial quantity of molecular gas persists in the outer 
disk, between radii of 40 and 200 AU, where photochem- 
istry from stellar and interstellar radiation dominates. 
The lack of molecular gas emission interior to this radius 
as indicated by our observations, combined with the lack 
of dust emissi o n wit hin a radius of 30 AU inferred by 
IWahhaj et all (|2007l ). implies that the 49 Ceti system 
appears to be clearing its gas and dust from the center 
out. The mechanism responsible for this central clearing 
is not indicated; in general, the best-developed theories 
to explain this transitional morphology are (1) central 
clearing through the influence of a massive planet and 
(2) photoevaporation by radiation from the central star. 

The clearing of gaps and inner holes has long been pre- 
dicted as a consequence of the fo rmation of massive plan- 
ets in circumstellar disks (e.g. iLin fc Papaloizoul 119861 : 
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iBrvden et alJ ll999). In the case of 49 Ceti, the forma- 
tion of a Jupiter-mass planet would be required at a dis- 
tance of ~ 40 AU from the star, roughly the inner ra- 
dius of the observed hole in the gas distribution. Such 
a scenario could also help to e xplain the size segrega- 
tion of dust grains observed bv lWahhaj et al.l (|2007f k a 
predicted consequence of inner disk clearing by gravita- 
tional influence of a massive planet is a filtration of dust 
grains by size, with only those below a certain threshold 
(typically 1-10 (im) accret ed across the gap along with 
a reduced amount of gas (|Rice et al.1 I2006D . However, 
this scenario ultimately requires the accretion of sub- 
stantial amounts of gas into the inner disk, a nd searches 
for molecular gas in the inn er disk of 49 Ceti (|Chen et al.l 
l2006t[Carmona et al.ll2007f ) have not detected such a pop- 
ulation. Another indication that an inner hole is likely 
caused by a massive planet in formation would be non- 
axisymmetric features resulting from its gravitational in- 
fluence, such as spiral waves. While the CO emission 
from 49 Ceti does not appear asymmetric within the 
limits of the SMA observations, more sensitive spatially 
resolved observations could address this hypothesis. 

The absence of gas in the inner disk is, however, consis- 
tent with a photoevaporation scenario: as the photoevap- 
orative wind produced by stellar radiation becomes com- 
parable to the accretion rate in the disk, material within 
the gravitational radius R g — GM±/ c 2 s will quickly drain 
onto the star, l eaving an evacuated inn e r hole free of gas 
and dust (e.g. iHollenbach et al.l 1 19941 : I Alexander et al.1 
l2006h . The gravitational radius for 49 Ceti is roughly 20 
AU, which is comparable to the inferred inner radius of 
40 AU for the outer disk. The larger outer radius may 
in fact be consistent with the later stages of photoevap- 
oration, after the inner disk has become optically thin to 
ultraviolet radiation and the inner disk radius slowly in- 
creases under the influe nc e of the photoevaporative wind 
([Alexander et al.ll2006[ K I Alexander fc Armitagd (|2007l ) 
propose a method of discriminating between inner holes 
caused by photoevaporation and those caused by the for- 
mation of a giant planet, involving a simple comparison 
between two observables: the disk mass and the accre- 
tion rate. As there is no measured accretion rate for 
49 Ceti, we cannot apply the criteria presented by these 
authors; however, we note that the low disk mass does 
indeed fall within the parameter space consistent with 
a photoevaporative scenario. Further observations are 
necessary to determine the origin of the inner hole; in 
particular, stringent limits on the accretion rate could 
suggest a photoevaporative mechanism. 

There are few disks which appear to be in a similar 
evolutionary stage to that of 49 Ceti; a rare example is 
the disk around the A star HD 141569. Like 49 Ceti, 
it hosts a disk composed primarily of sub/mi-size grains 
with infrared properties approaching those of a debris 
disk (jWahhai et al.ll2007t [Marsh et al.ll2002f) . while still 
retaining a substantial quantity of molecular gas with 
central region clear of gas emission, in this case out to 
a rad ius of ~ 11 AU (iGoto et al l 120061: IBrittain et all 
I2007t ). It exhibits a transitional SEP (|Merm et al.ll2004f h 
and observations of the rovibrational CO spectrum re- 



veal gas with disparate rotational and vi brational tem- 
perat ures (200 K and 5000 K respectively; IBrittain et al.1 
120071 ). indicative of UV fluorescence on the outer edges of 
an inner disk region cleared of gas and dust. An analysis 
of the chemistry and gas properties of the system similar 
to the one presented here for 49 Ceti was conducted by 
Uonkheid etafl (|2006l ). While the presence and extent of 
the inner hole are clearly indicated, the physical origin of 
this clearing is less obvious. The Br7 profile is ind icative 
of substantial accretion, and IBrittain et all (|2007f ) deem 
a photoevaporative clearing mechanism unlikely due to 
the large column density outside the cleared region and 
the la ck of evidence for a photo evaporative wind in the 
FUV (|Martin-Zaidi et all 120051) . However , iMerin et all 
(|2004h place a much lower limit of 10 -11 M /yr on the 
accretion rate, based on the assumed gas:dust ratio of 
100 and the low optical depth of the inner disk, which 
would be muc h more cons i stent with a photoevapora- 
tion scenario. IGoto et all (|2006l ) note that the rougn 
coincidence of the inner rim of the disk with the gravita- 
tional radius suggests that photoevaporation in concert 
with viscous accretion is a likely cause for the inner disk 
clearing. 

Whatever the origin of their morphology, the observed 
gas and dust properties indicate that the disks surround- 
ing both 49 Ceti and HD 141569 appear to be in a transi- 
tional state of evolution during which the dust properties 
are beginning to appear more like those of a debris disk, 
while the gas is in the process of being cleared from the 
disk from the center out. 

6. CONCLUSIONS 

The SMA CO J=2-l observations presented here pro- 
vide the first spatially resolved observations of molecular 
gas in the 49 Ceti system. The data reveal a surpris- 
ingly extended and complex molecular gas distribution 
in rotation about the central star, viewed approximately 
edge on and clear of molecular gas emission in the cen- 
tral region of the disk. Modeling the disk structure and 
chemistry in this system indicates that the inner disk 
is entirely devoid of molecular gas due to irradiation by 
the central star, while a ring of molecular gas persists 
between 40 and 200 AU, subject to photodissociation at 
the inner edge by stellar radiation. The disk model pre- 
sented here reproduces well the observed properties of 
the system, including the resolved CO J=2-l emission, 
the CO J=3-2 spectrum, and the spectral energy distri- 
bution. With dust properties similar to those of a debris 
disk and a substantial reservoir of gas maintained in the 
outer disk, 49 Ceti appears to be a rare example of a 
system in a late stage of transition between a gas-rich 
protoplanetary disk and a tenuous, gas-free debris disk. 
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